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Abstract: We show that the radion in a warped geometry bounded by two branes can 
have a potential suitable for inflation. Our construction is based upon a solution known 
in string theory as the linear dilaton, in which the back-reaction from a bulk scalar $ is 
exactly accounted for. The radion, stabilized by is much heavier than the TeV scale and 
its couplings to the standard model are much more suppressed than in the usual Randall- 
Sundrum solution. We present a new formalism for obtaining approximate time-dependent 
solutions, based on perturbing the exact solution to the coupled Einstein and scalar field 
equations in the bulk. It allows the radion potential to be computed directly in terms of the 
brane potentials for We show that simple exponential potentials on the branes can lead to 
a 4D radion potential with a flattened hilltop form, yielding inflation with a spectral index 
of typically = 0.96 and no higher than 0.99. With more complicated brane potentials, the 
descent from the hilltop can be a linear potential, giving a tensor-to-scalar ratio as large as 
r = 0.07 with Ug = 0.974. The couplings of the radion to the standard model particles are 
dictated by general covariance, so the details of reheating are explicitly calculable, leading to 
a reheat temperature of at least 10^ GeV. The quantum corrections to the inflaton potential 
from its couplings to matter are also calculable and are shown to be small, so that the 
prediction for the shape of the potential is under theoretical control, even with superPlanckian 
field excursions. 
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1. Introduction 

The possibility that we live on a brane embedded in a higher-dimensional space has been 
widely explored since one of its earliest string theory manifestations, the Horava-Witten model 
[P. Brane- worlds afford intriguing opportunities to go beyond the standard model of particle 
physics and cosmology, but they also present challenges for reproducing known physics. The 
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existence of moduli in these models is a generic feature with potentially important cosmo- 
logical consequences. If unstabilized, they can conflict with big bang nucleosynthesis or fifth 
force constraints. This was the case in the original proposal of Randall and Sundrum (RSI) 

which elegantly addressed the weak-scale hierarchy problem. These issues were resolved 
by introducing a bulk scalar field to stabilize the radion, the modulus determining the size of 
the extra dimension Q. 

Moduli are natural candidates to serve as the inflaton, and the radion is foremost among 
them. However it has been difficult to make radion inflation work in practice. Ref. |^ 
proposed the modulus of large extra dimensions |||] as the inflaton, but this required an 
inflaton potential with a rather peculiar shape to have inflation at the TeV scale Q. In the 
framework of warped extra dimensions, most of the effort has been either on inflation driven by 
motion of branes within the bulk, or else conventional inflation on one of the branes modified 
by an unconventional Friedmann equation due to the extra dimensions. We are not aware of 
previous proposals for getting inflation from the radion in the RS-I (orbifold bounded by two 
branes) scenario.^ It would be difficult to make this work, first because the Goldberger-Wise 
mechanism for stabilizing the radion does not give it a flat potential, and secondly because at 
Planck-scale values of the canonically normalized radion field, which could in principle have 
given chaotic inflation, the radion has already passed the barrier in the potential separating 
the compactified and decompactified regimes of the theory. In fact the light KK modes of 
the extra dimension become excited, invalidating the 4D effective description, long before the 
radion can reach the Planck scale. 

In this work we present a model of a warped 5D radion as inflaton that overcomes these 
difficulties. We incorporate a warped background solution that has negative curvature but is 
not AdS; instead it is conformally fiat, and relies on a linearly varying bulk scalar field for 
the source of its negative bulk energy density. The weak scale hierarchy problem is naturally 
solved, and inflation occurs at a high scale, with weak couplings of the inflaton to matter on 
the standard model brane. We are able to relate the potential of the radion to the bulk scalar 
field potentials on the branes in a rather simple way, and to find radion potentials (starting 
from brane potentials of exponential form) that are flat enough for inflation with a certain 
amount of tuning of parameters of the brane potentials. 

Generically, one always expects decompactified flat space to be a solution to the Einstein 
equations, and so the radion potential naturally has degenerate minima representing compact 
and noncompact fifth dimension, separated by a barrier. In our model the barrier is at 
Planckian field values, and we can find examples where superPlanckian excursions during 
inflation lead to an observable tensor component in the perturbations. Other examples have 
a small tensor component with inflation ending close to the hilltop. Because the minima are 
degenerate, we have the situation of topological inflation p, ^: the universe divides into 4D 
and 5D domains separated by eternally inflating domain walls where the radion is near the 
top of the potential. This guarantees the existence of regions where the inflaton is arbitrarily 

recent reference using different moduli fields is [^]. 
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close to the potential maximum (up to the limitation imposed by quantum fluctuations of the 
field) so that one does not have to rely upon fine-tuning of the initial position of the inflaton 
to get enough inflation. 

In section |2| we introduce the static solution with stabilized radion and contrast it to 
the RS solution, especially with regard to the hierarchy problem. In section |3| we perturb 
away from the static solution to find time-dependent ones, in the approximation that the 
back-reaction is small. We derive the effective 4D theory in section ^, computing the mass 
of the radion and the spectrum of Kaluza-Klein excitations for the radion and bulk scalar. 
We then construct some inflationary models using exponential brane potentials and work out 
their predictions for the CMB power spectrum, using both analytic and numerical techniques. 
Section ^ derives the couplings of the radion to the standard model, enabling us to compute 
the reheat temperature at the end of inflation and to estimate the radiative corrections to 
the inflaton potential. We give conclusions in section ^ and the bulk Einstein and scalar fleld 
equations in the appendix. 

2. Unperturbed static solution 

We begin by introducing the framework for our 5D brane-world model. The action for 5D 
gravity coupled to the stabilizing scalar field, $ = $ (y, t) with a bulk potential V (<I>) and 
potentials on the branes Vb,i (<!*) is 

-4 J d^x [i^]|,=o -4! Val [K]\y=^ (2.1) 

The extra dimension is an /Z2 orbifold (hence the factor of 2) and all functions are symmet- 
ric under y — >• —y. Thus we consider the interval y G [0, 1] for the compact extra dimension 
with the Planck (yo = 0) and TeV {yi = 1) branes located at the orbifold fixed points. We 
also include extrinsic curvature terms [K] that are needed for properly defining the variations 
of the action. 

The most general metric ansatz that respects 3D homogeneity and isotropy can be written 

as 

ds' = e^^^^'yUe - e2^(*'?^)dx2 - e2^(*'^)dy2 (2.2) 
The field equations for this ansatz are well known and we reproduce them in appendix 
2.1 Static solutions 

It will be useful for what follows to construct an exact solution to the coupled bulk scalar 
field and Einstein equations, which is conformally flat: 

ds'' = ^<y^dt^-d^^-b'dy^) (2.3) 

Notice that this differs from AdS5 in Randall-Sundrum coordinates where the e^"^^) factor 
does not multiply dy^. The negative vacuum energy density provided by the bulk scalar 
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potential is not constant, so the geometry is not AdS. This solution is well-known (see e.g., 
[^]-|12]) in the context of strings propagating in a spacetime with subcritical dimensionality, 



compensated by a spatially linearly varying dilaton similar to our bulk scalar. With this 
ansatz the Einstein equations in the bulk reduce to 

-3 {n" + = Qci>'2 + 52g2n y ^^^^ (2.4) 

-6 (n) ' = 4 (^-^$'2 + ft^e^" V (c^)) (2.5) 
As usual, the scalar field equation of motion 

$" + 3nV-62e2"^ = (2.6) 
is not independent of the Einstein equations, but can be obtained by taking linear combina- 



tions of (2.4, |2.5| ) and derivatives thereof. By choosing the bulk potential to take a special 
form, V = — y exp(cy$), one finds that the solutions for n and $ are simple linear functions 
of y: n{y) = no — khy, ^{y) = + c^y. Here k plays the role of the AdS curvature scale, 
as in the usual RS solution, and by analogy we will so refer to it, even though the solution 
is not AdS. The integration constants no and $0 can be absorbed into the normalizations 
of b and V respectively without loss of generality, so we set them to zero. The equations of 
motion then imply th 
potential is given by 



motion then imply that V = — cy = 2-^, = With these choices, the bulk scalar 



where we define 



and the solution is 



V{<^) = -^l,lk^exp{2f,-'<^) (2.7) 



\/3 , , 

fi5 = — (2.8) 
H5 



n{y) = -kby; $ = fi^kby (2. 



In contrast to the original RS construction, we do not include a bulk cosmological constant; 
the negative bulk potential V{^) is responsible for the curvature of the 5D geometry. Similar 
negative potentials have been studied in the context of Golberger-Wise stabilization of bulk 
fields in AdSs [^], where it was observed that the unboundedness of the potential does not 
lead to instabilities in the AdS background The present case is similar; in the analogous 
subcritical string theory situation where the bulk scalar is the dilaton, in the string frame the 
negative coefficient of the potential is interpreted as a negative bulk cosmological constant, 
which would give rise to an AdS background if the dilaton were held fixed. In the present 
situation we have in addition the couplings of the bulk field to the branes, which prevent the 
field from running away. 
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2.2 Junction conditions 

The Israel junction conditions (see e.g., ref. |T^) and the boundary conditions for the scalar 
field are given by 



^Vr-^ 3 



(2.10) 



Vi 



where the ± apply respectively at yo = and yi = 1. Then explicitly, 

2^4k = Fo(0), 2iilk = - Vii^i) (2.11) 

2/.5fc = ||(0), 2^5fe = -e"^||($i) (2.12) 

where e"'^ = e~^^ and $i = pL^kb are respectively the warp factor and the scalar field value 
on the TeV brane. 

Let us recall the physical significance of the junction conditions. In the pure RS model 
with no scalar field, the two junction conditions require fine-tuning of the brane tensions (here 
represented by the values of the potentials Vq and Vi on the branes) . For generic values of the 
tensions, one would obtain a bulk solution that is not static, and which contains a black hole 



in the extra dimension |16]. Thus the two tunings can be interpreted as (1) the usual setting 
of the 4D cosmological constant to zero, to obtain a static solution, and (2) the tuning of the 
bulk black hole mass to zero. The radion h is exactly massless in this solution, and can thus 
take any value. 

Next consider the addition of the bulk scalar field. In the usual implementations of the 
Goldberger-Wise mechanism, the back-reaction of the scalar field on the metric is taken to 
be small, and one simply solves the bulk scalar equation in the AdSs background. This is a 
second order equation, and the two constants of integration in the solution are determined by 
the two additional boundary conditions involving dVi/d^. This generically induces a potential 
for the radion that is minimized at some value 6o- 

In our solution, the situation is somewhat different; we have singled out a particularly 
simple form of the scalar field solution, which is only compatible with certain choices of 
dVi/d^^ at the boundaries. It is still true that the radion is stabilized as in the generic GW 
mechanism, but there is an additional tuning of potential parameters needed to maintain the 
linear solution for This extra tuning is on the same footing as that of the bulk black 
hole mass in the pure RS solution. We suspect that our results can be generalized to more 
complicated solutions in which this tuning is relaxed, but for ease of computation, we will 
adhere to this special situation. 

The upshot is that the junction conditions amount to three tunings of brane potential 
parameters, plus one condition that fixes the value of the radion h, and hence the warp factor 
on the TeV brane. Let us illustrate with a simple example of brane potentials which is similar 
to that made by GW: 

yo = ?T^g($ + no)^ Vi = -m2($ + ui)2 (2.13) 
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Solving the junction conditions, we obtain three constraints that can be regarded as fine- 
tunings on the parameters tuq, uq and ui, 

k 

1^1 = 2 ' '"0 = = (2 - kb) (2.14) 

while mi can be regarded as adjustable, and its value determines the warp factor through 

e-'^"" = ^ (2.15) 
nil 

One notices a disadvantage of our approach relative to the usual one; we have had to build 
the hierarchy of scales into our original choice of brane couplings to get an exponentially 
small warp factor, as opposed to having generic Planck-scale values for the brane potential 
parameters. The next example shows that this problem can be ameliorated if the potentials 
take an exponential form. 

2.3 Exponential brane potentials 

In our subsequent construction of an inflationary solution, we will make use of a different 
choice of brane potentials, which are exponentials in 

Vo = Ao e°o*/'^5 - Ao e^o*/'^^ (2.16) 
Vi = -Ai e°i*/'^5 + Ai e''i*/'^5 (2.17) 

The junction conditions lead to 

Ao = 2^4k , Ao = 2,lk (2.18) 



Ai = 2/x^A; e^^-^i^'^"", Ai = l^k { \ e^'-^'^^'^ (2.19) 

V«i-Pi/ \ai-pij 

For a given desired value of bo determining the hierarchy of scales between the two branes, 
there are thus four free parameters Oj, /3j that can eventually be used to tune the potential 
of the radion to a form suitable for inflation, while the dimensionful couplings Aj and A, are 
fixed in terms of these. Notice that if ai and /3i happen to be moderately close to 1 {e.g., 
~ 1.1), the explicit hierarchy can be ameliorated between the dimensionful parameters on the 
Planck brane versus those on the TeV, while yielding a sufficiently large value of kbo ~ 37 to 
explain the TeV scale, since then Ai/Aq ~ e(^-"i)''^o and Ai/Aq ~ e^^^^i^'^^'o. 

3. Time-dependent solutions 

Our goal now is to extend the static solution of the previous section to dynamical ones in 
which the radion is displaced from its stable equilibrium value 6o- Ii^ refs. [15|,[^]-[20], 5D 



solutions that also include the dynamics of a bulk scalar field were found for the case of a 
single brane. In these papers the search for solutions was simplified by assuming that the 
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metric functions and the bulk scalar were separable functions of t and y, which we will justify. 
The key to our approach will be to allow for a large excursion of the radion field (which plays 
the role of the inflaton) in response to linear perturbations of the metric and bulk scalar. 
Although such an ansatz is not guaranteed to be consistent, we will show that for suitable 
choices of the parameters for the exponential brane potentials (2.16,2.17) it is justified, and 
will allow us to find physically interesting inflationary solutions. 

3.1 Small perturbations with large radion fluctuations 



We make an ansatz for general scalar perturbations similar to that in ref. |21], 



$(x,y) 



in{y)+ip{t)) + 6^x,y) 



(3.1) 
(3.2) 



It reduces to the static solution of the previous section when F = G = 6^ = ip{t) = 0. 
Normally then, one would linearize in all of these quantities. However we will find a more 
general solution for the perturbations by treating only F, G, 5^ as being small, while working 
to all orders in {p(t) (and (p). This corresponds to allowing for large excursions of the radion 
during inflation. We will justify it a postieriori by choosing special values of the brane 
potentials that lead to a sufficiently flat inflaton potential. 

The perturbations F, G, 6^ are constrained by the nondynamical Einstein equations, 
namely those containing only first derivatives. The G05 equation gives 



(G + /U5 M$)n' + {F' + fi^^6^')^ -F' = 
Similarly the off-diagonal Gij equation gives 

2didjF + didjG = 

It can be integrated to find 

G{x,y) = -2F{x,y) + Go{t,y) 
where an untenable term linear in x has been set to zero. The Gqj equation 

{djF + iJ^^djS<^)^ = 



can be solved for the bulk scalar perturbations, using eq. (3.5) : 

6<^{x, y) = -H5F{x, y) + (5$o(i, v) 



Lastly the Gjs equation, after substituting eqs. ( |3.5| ) and (3/?), reduces to 

dj{F' + ?>n'F) = 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 
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The general solution to ( |3.8| ) is 

F{x, y) = F(x)e-3"fe) + Fo(t, y) (3-9) 

We find that the F{x) terms in the solutions to the metric and scalar field perturbations, eqs. 
( |3.5| , |3.7| , |3.9| ), vanish at linear order from the off-diagonal constraint equations and so F{x) 
is undetermined. Therefore it can be consistently set to zero. 

Now let us specialize the general ansatz to a form that is adapted for the cosmological 
solutions we seek: 

=n(y)/Mt)), F^i^^y) = h{^{t))e^^+^)^^--<y)) (3.11) 
$(t, y) = /X5 {kh^y - v{t)) + 5^{^{t),y) (3.12) 

where ni = n{l) = —kbo and e is an adjustable parameter.^ It reduces to the previous 
static solution when the radion fluctuation (p(t) vanishes, if /, /2 and 6^ also vanish. One 
can interpret /, /2 and 6^ as small back-reactions induced on the metric and scalar field, 
respectively, by the (possibly large) radion fiuctuation ip. This ansatz does not provide an 
exact solution to the equations of motion, but it is a good approximation at leading order 
in /, /2 and 5$ as long as these quantities are small, in a sense to be specified. We will 
find situations where it is possible to reliably consider large excursions of the radion during 
infiation, even superPlanckian ones. It will be shown that the functions /, /2 and 6^ can 
be algebraically derived from the brane potentials Vq and Vi. In turn they determine the 
potential for the radion in the 4D effective theory. We will find that the radion potential 
derived from this ansatz is dominated by contributions from / and its bulk scalar counterpart 
(the part that is fiat in the extra dimension). The other parts of the perturbations, those 
with exponential y-dependence, turn out to be relevant solely for the stabilization of the static 
solution. 

3.2 The Go5 equation and junction conditions 

To arrive at the functional form of the scalar perturbations we examine the off-diagonal 
Einstein equation. Linearizing the left hand side we find 

_^(l + /_|_3F. + ^)=,-«' (3.13) 

This motivates us to parametrize 6^{ip,y) = + 5^2 where 

f,^'6<f,i^,y)=gi^), ^5ij<I>2(v.,2/) =52((^)e(='+^)("i-"(j'» (3.14) 

This separation is useful, because it will be necessary to require that /, /2,9,52 1 but not 
necessarily that 93 <C 1. With this parametrization, the dominant "1" terms which are zeroth 



^We will see that the choice e — —1 leads to small back-reaction of /, /2 and 5$ while allowing for large ip 
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order in the perturbations cancel as in our previous analysis and the linearized Go5 equation 
becomes 



It is satisfied everywhere in the bulk if 

dip dip 



f 



, .df2 dg2 /o , , \ 

(l+e)^ + ^ = (3 + .)(/2 + <72) 



and has negligible y-dependence when 



and /2,52<1 



(3.16) 
(3.17) 

(3.18) 



In our numerical analysis we find eqs. ( 3.16| , p. 17 ) hold to a high degree of accuracy when the 
back-reaction is small. Moreover, for the cases we study, eq. ( 3.17] ) will reduce to eq. ( p. 16 ). 
Now let us reconsider the junction conditions in the presence of the perturbations: 



24k 1 + / - (1 + e)/2e- 



-(34-£){ni-ni) 



2^^k (l + ((3 + e)(72 + 2/2)e-(3+^)("»-"i)) e""' 



y=yi 



y=yi 



y=yi 



y=yi 



(3.19) 
(3.20) 



where now = —/is 99 and <I>i = n^ikbo — take the place of the unperturbed scalar field at 
y = and y = 1. These are algebraic equations that determine the back-reaction functions 
/, g and /2, 52- At low energies the functions all have a linear dependence on ip and are related 
as g = f and g2 = (1 + e)/2. However, when e = —1, g2 vanishes at linear order in ip so 
that the bulk scalar perturbations are effectively independent of y. Even when ip is large, this 
choice simplifies the subsequent analysis without affecting the radion dynamics, so we will 
adopt e = — 1 from now on. 

The linearized jump conditions can be solved to first order in g by Taylor-expanding Vq 
about $ = —fJ.5<p at y = 0: 

2ke-^ (1 + (2/2 + 252)e2"i) - f,^^ 



a 



-526^"^ + 



V" 



(3.21) 



where the prime on potentials denotes Expanding next for Vq in eq. ( 3.19| ) we find 



/ 



-1 + 



2/iiA; 



(Vb + M5(9 + 52e2"^)yo') 



(3.22) 



Similarly, at y = 1, we expand V{ about ^ = <I>i — fi^ip and find for arbitrary stabilizing 
potentials using (|]2l|) and ( p2D 



h 



-92 



1 Vo"yi"e"ie'^(l - e2"i; 
2^ 4A;(yo" + y|"e3"0 



(3.23) 
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Then expanding for Vi in eq. ( 3.19 ) 



where 



G2 



(3.24) 



(3.25) 



Vie is tempting for simphfying 



A fine-tuned choice of brane potentials such that Vq = 
the solutions since then 52 = /2 = 0. However, our study of the KK excitations (section 4.2) 



shows that this choice is precluded by the presence of a tachyon instability. The expressions 
for the back-reaction functions may nevertheless be simplified by neglecting the subdominant 
contributions involving extra powers of the warp factor: 



9 



f 



92 



f2 



-1 + 



V" 



(^0 + /i55^o) 



2A:e-'^^5(l + /) + /^5'^ie"^ 



Vie 



-92 



1 + 



1 



^-7r(/^5-Me-+<e-,) 



(3.26) 
(3.27) 

(3.28) 
(3.29) 



2 4A; 

Our results for f,g and /2,S'2 are derived from the jump conditions and not from the 



Einstein equations, so it is not obvious that they satisfy the relation ( 3.15 ) from the G05 
equation. However, it is simple to show that this indeed is the case, at linear order. Differ- 
entiating the first jump condition (at y = or y = 1) in ( 3.19 ) with respect to if using ( p. 14 ) 



gives 



dV 



y=yi 



dip dVi 

dip 



^5 



y=yi 



dip dip 



dip 



(3.30) 



y=yi 



Eliminating ^ using the second jump condition (at y = or y = 1) in ( p.2Cl| ) and linearizing, 
the dominant "1" terms again balance and we arrive directly at ( 3.15| ) evaluated on the 
branes (which also guarantees that ( |3.15 ) is satisfied in the bulk since both sides of the 
equation vanish). This is similar to what was found in ref. ||2^ where the linearized off- 
diagonal Einstein equation, G^s evaluated on the branes was shown to be equivalent to one 
of the junction conditions of the metric so that no new constraints arise. 

For small 99 ^ 1 we can Taylor-expand about the background and arrive at simple 
expressions for the back-reaction functions 

-gp {ip - gae-'^''") apip 
9 = 1, 9 = (3.31) 

1 - CTo 

o-o) 



1 



92 



0, /2 



ip{ai 



1-ao 



-2fci,o(l_^^)) 



1 y(tTl 

2 



1-ao) 



(3.32) 
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where the jump conditions 



±/i5-ie"-l^/(^.)|$, = ±e^^^ifV,m^^ = 2k (3.33) 

have motivated us to parametrize the brane potentials such that 

±e"»y/'($)|$, =2fc(l-a,) (3.34) 

This has a convenient physical interpretation. When cji = the radion is massless and the 
background solution will be BPS in the sense of the solution generating technique utilized in 



ref. |13]. In general, the depend on parameters of the brane potentials, and they determine 
the zero-mode masses for the radion and the bulk scalar field in terms of the curvature scale 
k. The effects of various choices of fij will be explored in section 

3.3 The remaining Einstein equations 

So far we have succeeded in finding the back-reaction as a function of the radion fluctuation 
99, but we have not yet determined the dynamics of <p. This comes from solving the remaining 
Einstein equations. They are difficult to solve exactly due to their explicit dependence on 
y through the warp factor n{y) which also appears in the metric and scalar perturbations. 
However, when we supplement the restrictions (|3.18 ) on /2 and 52 by analogous ones for /, 



f,g<.ip and f,g<.l (3.35) 

then it is possible to linearize to obtain some equations that have negligible y-dependence. 
The radion will induce only a small back-reaction on the 4D slices of the 5D spacetime relative 
to the static solution, 

dsl ~ e2"(f) {dt^ - a^d^^) (1 + 2nf (99) + 2/2(vp)e-2("-"i) + . . .\ (3.36) 



Similar to the off-diagonal equation ( |3.15| ), zeroth order terms vanish since these involve 
factoring out 99 dependence with the background solutions. The remaining equations can 
then be averaged over the extra dimension to find the radion dynamics: 

/' dy^G^ = 4 [' dy^T^ (3.37) 
Jo Jo 



It was shown in reference [23| that the cosmology of the radion in the RS model can be 
obtained equivalently from the effective action or from averaging the Einstein equations as in 
( |337| ). In Section 4.1 we will construct the 4D effective action that leads to the bulk-averaged 
equations of motion. 



The bulk equations to be averaged are given in the Appendix, eqs. ( A.2D -( [A~^ ). Collecting 



t derivatives to the right hand side, Using eq. ( |2.7| ) for V and ( 3.12 , 3.141) for linearizing 
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eq. ( 3.37 ) with Kg 
functions 



3^5 and integrating over y we find to first order in the back-reaction 

2 



- +-(^ 1 + Ti 
a J a 



a 

a la 
-+ - 
a \a 



.a 



2- + ( - ) + ( 2-(^ + ) (1 + Ti) + ^ + T 



(3.38) 
(3.39) 



+ — 1 + ^ + 2T4 

2 I 2 



A:2e-2<P (4/ - 35 - Te) (3.40) 



where primes denote ^ 
2 



and 



Ti = 



T2 EE --Ti - 5' - 36-3^"° (/^ + g'^) (1 - e'=^«)f) (3.41) 
T3 = -T2 - 45' - 12e-3'=''« (/^ + g'2) (1 - e'^^")^^, = Se'^''^'' f!^{l - e'=^°)f] (3.42) 



526 



'3kbo I 



(1 - e^^°)J7, Te = 3e-^'=''o (2/2 + 552)(1 - e'''°)n (3.43) 



with 



Jl = (l-e-3^'^n)"^ (3-44) 
All of the Tj terms are subdominant to other similar terms appearing in the equations, so 
it is consistent to drop them and keep only the explicit / and g source terms (notice that 
Ts^e are exponentially suppressed by powers of the warp factor). This demonstrates our 
earlier assertion that the radion dynamics are dominated by the Planck brane potential, 
which determines / and g, while /2 and g2 are needed for stabilization only. The remaining 
Einstein equations then take the simple FRW-like form 



a\ a . 
a J a 

a 



^ = 3A:2e-2^ (/ - g) 



2-+ - +(p + 2-^ + -^'^ = 9k^e-'^'P{f -g) 



a- 2 

2 .;,2 



^+(-) +^ = ^V^-(4/-3,) 



These equations are augmented by the linearized Klein-Gordon equation 

(p + 2,% + ^^ = Ak^e-'^ (f-\g 



(3.45) 
(3.46) 
(3.47) 

(3.48) 



found at this level of approximation by integrating (A.l) or equivalently by starting from the 
4D effective action (4.3) in Jordan frame which we derive in the next section. 

Of course, not all four of the equations ( |3.45| - p^48| ) are independent. As was observed 
in ref . |Q , the linear combination of equations ( |3.48| ) — ( |3.46| ) -|- 2 x ( |3.47| ) is equivalent to 
eq. ( lOSi ). Furthermore, f (^^+^^(3 x ( |3l5D - <^M^)/v + ( p8| ) + (^45]) is equivalent 
to eq. ( ^.471) . The two constraints arise as a consequence of the gauge symmetries, namely 
reparametrizations of the t and y coordinates |22|,[23|. 
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4. Inflation in the effective 4D theory 



It is not difficult to derive the form of the effective action that gives rise to the 4D equations 
of motion ( |3.45| )-( |3.48| ). Starting from the 5D action (^]^) and the metric ( |2.2| ) one arrives at 

"1 



-J d^Xy/glVo\y^Q- j d^x^Vily^^ 



(4.1) 



The boundary contributions from the compact extra dimension exactly cancel the extrinsic 
curvature terms and so do not appear in this expression. Substituting for the metric functions 
A, B and N in (|3lo| ), for V and $ with and (Isl^ ), and using the jump condition ( ^.191) 
to rewrite Vi in terms of the perturbations one can linearize about <I> = ^^{kboy — ip). After 
integrating over y we arrive at the effective action 



2/x: 



3kQ 



a^e'^dt 



-j (1 + Ti) + ^2 Q + T2 ) - 3fc2e-2'^ U-g- T„ 



(4.2) 



Neglecting the small Tj terms as before, we can write the simplified effective action in Jordan 
frame as 

1 



S 



2k^ 



(4.3) 



where V, 



r,J 



'•2^4 ^"^(Z — g) and the 4D Newton's constant K4 = SvrG is found from 



K4 — 



-3kboy 



dy 



9kQ 



(4.4) 



It is straightforward to show that this action implies the equations of motion ( 3.45| )-( |3.48| ) to 
linear order in /, g when we make use of eq. (|3l^ ). 

One can go to the Einstein frame by performing the Weyl transformation 



a (i) = 7a (r) , dt = -)dT, 7 = 6 2 
after which the action takes the form 



5 



2kI 



dr a-^ 



6 - + - 2KiVr{ip) 



(4.5) 



(4.6) 



where dots now denote and the radion no longer mixes with the scale factor. The Einstein 
frame potential is 



Vr 



9k' 



if -9) 



(4.7) 



and (f is related to the canonically normalized radion field <j) by 



(4.8) 
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where mp is the reduced Planck mass. The Friedmann equation then takes the usual form 



3 V2 



+ Vr 



as does the Klein-Gordon equation 



(4.9) 



(4.10) 



We will continue to do some of the subsequent analysis using the dimensionless field (p however, 
to avoid having to repeatedly write ^4,. 

4.1 Radion stability 

To summarize the results to this point, we have shown that in the class of 5D warped models 
we are considering, the radion potential is given by 

VM) 



if -9) 



(4.11) 



where the functions / and g are determined by the Planck brane potential Vq to leading order 
for small /, 5, by eqs. ( 3.26 , 3.27|) . The TeV brane potential Vi has virtually no effect because 
of the suppression of /2 and g2 by powers of the warp factor in eqs. ( 3.21| - 33^ ). 

We are assuming that = is a stable equilibrium point. Since f = g = f — g' = at 
(j) = (see eq. ( 3.16| )), it is clear that = is a critical point of the potential, and Vr also 
vanishes at this point, leading to a Minkowski solution. To compute the radion mass, notice 
that /" - g" = f. Therefore 



0=0 



2fcV'(0) = 2A;V(0) 



(4.12) 



% and ^o" 



d¥2 



Using eq. ( ^) , we find that ^'(0) = 1 - 2k/V^'{Q). (Recah that g' 

Consider for example the exponential potentials ( 2.16| - ^T7 ). The radion mass is deter- 
mined by V^{<d) = 2k[al{l - /3o) - - ao)]/(ao - /3o)- Recalling our discussion of the 
hierarchy problem, where it was pointed out that values of ai and /3i close to unity alleviate 
the need for a strong explicit hierarchy between the scales of the potentials on the two branes, 
we are motivated to define 

ai = l + a, /3i = l + /3 (4.13) 



and for consistency when comparing parameters 



/3o = 1 + /3 



(4.14) 



Then it is straightforward to show that ^'(0) = q/3(1 + a/3)/(l — a/3). The radion mass 
squared is thus given by 



-2k^aP 



1 + a/3 
1 - a/3 



(4.15) 
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and there is a second zero mode whose mass has a similar expression in terms of the a and 
/3 parameters. Hence is generically nonzero, and it is positive as long as — 1 < a(3 < or 
a/3 > 1. As for the remaining parameters we will show in the next subsection that < af3 <^ 1 
is required in order to avoid a tachyon instability in the second light state, and to address 
the hierarchy problem. 

We emphasize an important difference between our model and the RS model supple- 
mented by Goldberger-Wise stabilization of the radion. In the latter, the radion mass is 
suppressed by the warp factor and so is naturally at the TeV scale, whereas in ours, whose 
background geometry is different, it is natually at the scale k; we must do one tuning of 
parameters, a <C 1, to make it small enough for inflation. This will be advantageous when we 
consider the robustness of the radion-as-inflaton potential against quantum corrections from 
standard model physics later on. 

4.2 Kaluza-Klein excitations 

One might wonder if it is sufficient to stabilize the radion without also considering the zero 



mode of the bulk scalar field. As ref. |21] has shown, there is only a single KK tower for the 
two fields. However, in our model there is a second zero mode in the tower, the sign of whose 
mass squared depends on the choice of brane potentials. We will show that this state can 
be made much heavier than the radion so that it is consistent to neglect it along with the 
higher KK modes. Our main point is to show that the nonzero KK modes are not excited 
during inflation, which if it happened would invalidate our effective 4D description. We will 
find that these modes have masses of order k, while the second zero mode hass mass -Jaik, 



with (Ji defined in ( 3.34 ). The details of the calculations are left to Appendix^. 



Following Csaki, Graesser and Kribs (CGK) [^T| we consider the spectrum of pertur- 
bations around the background solution. The metric ansatz used to describe the scalar 
fluctuations differs slightly from theirs since we are working with a y-coordinate in which the 
background is conformally flat 

ds^ = e^<y) [e^^^^^'y^ri^^dx^dx^ - e^^^^'^Uy^] (4.16) 
^(x, y) = $6(y) + 5^{x, y) = fi^kby + 6<^>{x, y) (4.17) 

A similar analysis was subsequently done by Kofman, Martin and Peloso (KMP) [^]. They 
generalized the CGK results by constructing a gauge invariant combination of scalar per- 
turbations which ensures the hermiticity and orthogonality of modes in the mass spectrum, 
defined here as 

v = z[F+ — (4.18) 



/ 

with 

z = V2-^e- = -V2/i5e-— (4.19) 



n' 

V is analogous to the Mukhanov- Sasaki variable ^ in 4D inflationary cosmology. Ex- 
panding the action (|2.1| ) to second order and diagonalizing in terms of the different scalars. 
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one finds separable solutions v{x,y) = J2i Qj{^)^j{y) such that the quadratic action of the 
KK modes can be reduced to 



with the normalization coefficients 



S = ^Cj j d^x Qj [-□ - m] ] Qj (4.20) 



2/0=0 



where Fj is related to vj as F is to in (4.18). The solutions are determined by the jump 



conditions ( |B.6| , |B.7| ) and the two nondynamical constraints ( [B.3| , |B.4| ) originating from the 
linearized off-diagonal Einstein equations to determine G and 6^. 

To calculate the mass spectrum in our model we incorporate the results of KMP and 
CGK except in the application of the boundary conditions. Unlike the present work, refs. 
[HUH] assumed a convenient stiff potential limit for the stabilizing potentials that forces 
6^ to vanish on the branes with the result that only the radion zero mode appears with the 
infinite tower of KK excitations. Relaxing this assumption gives rise to the bulk scalar zero 
mode. One can also understand the two light modes as being the moduli of the positions of 
the two branes in the bulk [p7| . 

General solutions obtained from the action expanded at second order for vj are equivalent 
to those obtained for Fj from the linearized Einstein equations [^]. We prefer to solve for 
Fj similarly to CGK since the boundary conditions and normalization look simpler in this 
language. The heavier KK excitations and the zero modes take the form 



where 



Fj = e-^^^'^y [Aj sini^jboy) + Bj cos{y/Xjboy)j (4.22) 
F, = e-i^^'^y i^A-.e^^^'^y + B-^e'^^^'^y'^ (4.23) 



A,=m|-^ , A. = ^-m^ (4.24) 



for 4m| > 9fc^ and ^ A;^ respectively. The index j > is an integer which labels 
the heavier KK modes, while z = r,s stands for the two light (radion or bulk scalar) zero- 
modes. The complete solutions are obtained by applying the linearized boundary conditions 
to determine the ratio of integration constants and the mass. 

The lightest KK excited state has Aq = with mg = |/c^. Heavier modes must be 
found numerically except for some special cases which we display in Table |^. Recall from the 
discussion in the previous subsection that ai = aP ^ 1 is needed to explain the hierarchy 
between the Planck and TeV branes. If |(To| <^ \cri\ <^ 1 so as to achieve a large mass gap 
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rtij (KK modes) 


(zero modes) 


(To < , (Tl > , 
|cro| < |c7i| < 1 


ml = |A;2 


U + h) 


2 


m^g ~ — 2croA;^, diA;^ 













Wi\ > 1 








2A;2 


ctq = (7i, \cri\ 1 











Table 1: A summary of the mass spectrum for various brane potentials parametrized 

by ai. The entries correspond to the general case where tachyons are absent, followed by the massless 
radion, the stiff potential limit and fine-tuned brane potentials. 



between the radion and the other modes, then the KK masses are well-approximated by 

1 2 



m^.+i = -k +k 



l\ vr 



^ 2 j khQ 



(4.25) 



for j > 0. Thus all KK mode masses are at least order k, independent of the precise details 
of the stabilizing potentials. 

We previously showed that the radion mass depends on Vq, the potential on the Planck 
brane; see eq. ( 4.12D ). Similarly the bulk scalar zero mode mass depends on the TeV brane 
potential Vi. In the case where the brane potentials are tuned such that \ai\ <^ 1 (eq. (|3.34D ) 
the zero mode masses are approximately 



m^ 



-2o"oA;^ , aik^ 



(4.26) 



To avoid tachyonic instabilities we must impose ctq < and ai > 0. This result agrees with 
the analysis of criteria regarding the presence of tachyons for brane world models utilizing 



the Goldberger-Wise mechanism carried out in reference [28|. 

Comparing ( [4.26 ) with ( [4.15| ), for the exponential brane potential Vq ( p. 16 ) which implies 
(To = we corroborate at first order in iTq the more exact expression for the radion mass 
given in the previous subsection. We will henceforth assume that ai = a/3 from Vi ( p. 17 ) is 
large enough so that the bulk scalar zero mode is much heavier than the radion: 



|(To| < (Ji « 1 



(4.27) 



We will verify that ^ H (the Hubble scale during inflation) so that it is consistent to 
ignore the bulk scalar zero mode during inflation. The KK modes are much heavier and thus 
can also be considered as frozen. 
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4.3 Radion inflation: an explicit model 



Let us now consider the full radion potential (valid to leading order in /, g) corresponding to 
the brane potentials ( p.16 - 2.17 ), with the restrictions ( 2.18 - 2?ig| ) and the definitions ( [4.14D . 
Ignoring the small exponential terms involving /2 and §2 we find that 



9 



(g - /?)e^^ + /3(1 + a)e(^-")^ - a(l + /?) 
a(l + /3)2 - /3(1 + a)2e(/3-a)'/' 



(4.28) 



and the combination (/ — g) that enters the radion potential ( 4.11 

1+5 



IS 



f-9 



a — /3 
9 



a — j3 



/? - a + a(l + /3)e" 



+ a)e 



-aip 



(4.29) 



In the case we are interested in where |a|, ^ 1 in ( [4.13| ), the expressions for /2 and 52 
( |3.28 , 3.29] ) take the simple form then 



/2 ^ + 9 + ^a/3((/9 - vj2) 
92 ^ f - 9 + ^aPif'^ 



(4.30) 
(4.31) 



Since ~ 1 during inflation while d/3 <C 1, it follows that if / and g are small then so are /2 



and g2 and the consistency conditions ( 3.18| ) are satisfied. 

To design a suitable potential for inflation, we have the two dimensionless parameters 
a, (3 and one dimensionful scale k at our disposal, subject to the restriction ( [4.27| ) on o"o = a(3. 
In addition ui = a/3 must satisfy ( 4.27 ) and we must ensure that f,g <^ 1 for consistency 
of the effective theory. Notice from ( 4.28| ) that if \a\ ~ |/?|, then generically g ~ 1 even if 
|a|, <C 1. However if 

\a\ < (4.32) 



(or f3 <^ a, giving qualitatively similar results) then / and g can be made parametrically 
small: both vanish in the limit a — t- 0. This motivates our choice of parameter values. We 
have also remarked that the radion potential has a second degenerate (Minkowski) minimum 
because of the decompactification limit as 99 — )• 00. Thus the generic form of the potential 
is that of a hilltop model with a barrier separating the minima at 99 = and 99 = 00. We 
find that this barrier can be made very flat when /3 is close to —3, for < a ^ 1, making 
it possible to get inflation from near the top of the barrier. The shape of the potential is 
illustrated for a = 10^^ and a few values of /3 near —3 in figure |l[ Figure ^ shows that 
the condition f,g<^lis satisfied in the region < 99 < 3 where inflation starting from the 
top of the barrier might take place. Choosing a = 10^^ and /3 = 6 x 10^'^ for the fiducial 
values of a = 10~^ and /3 = —3 is sufficient to ensure eqs. ( [1.30 , 4.31 ) for /2 and g2 are 
good approximations and that the mass squared of the second light zero mode is much larger 
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Figure 1: Radion potential for a = 10 ^ and 
P = -2.9, -3 -3.1. 
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Figure 3: 1^9 as a function of number of e- 
foldings for a = 10~^, (3 — —3 potential, 
starting Aip — 0.01 from the top of the barrier. 
Inset shows end of inflation. 



Figure 2: Logarithm of the functions g and 
f—g that must be ^ 1 in the region of validity 
of the potential, for the case a = 10^^ and 
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Figure 4: Spectral index as a function of 
number of e-foldings until the end of inflation, 
corresponding to fig. 3. 



than that of the radion. We will presently show that it also exceeds the Hubble scale during 
inflation. 

We have numerically integrated the inflaton equations of motion starting close to the top 
of the barrier. Taking the number of e-foldings N as the time variable, these can be written 
in the first-order form, 

dN ix^H' dN ^i^H dip ^ ' ' 

where it = (p\s the canonical momentum and the Hubble parameter is given by = ^vr^+IV^ 
in Planck units. (Recall that /X4 arises in the relation between the canonically normalized 
radion (j) and the dimensionless one 99. ) For a = 10"^, /3 = -3, the maximum of the potential 
is at ifm = 3.01. Starting at = ipm — Sip with Sp = 0.05, for example, gives Ng = 685 
e-foldings of inflation, p as a function of is shown in figure ^. Although the total number 
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of e-foldings depends upon how small is taken to be, the spectral index at the time of 
horizon crossing (nominally 60 e-foldings before the end of inflation) is insensitive to the 
initial condition, as long as at least 60 e-foldings of inflation occurred. This behavior has 



been observed in a number of previous studies of hilltop inflation models |29|-[p4|- 

To compute the power spectrum and spectral index we have used the prescription (see 
eq. (43) of ref. [||) 

in terms of the slow roll parameter e = ^mp(l^'/yr)^ = These are to be 

evaluated at horizon crossing, approximately 60 e-foldings before the end of inflation. More 
precisely, we take the crossing of the scale relevant for the COBE normalization of the power 



spectrum to be A^cobe = 56.5 [31|, a number that depends only logarithmically on the scale 
of inflation; since the spectrum changes slowly with N , this is an adequate approximation for 
our purposes. 

Fig. ^ shows the spectral index Tig 3jS 3; function of Nf> — N, the number of e-foldings until 
the end of inflation. We see that Ug = 0.96 at the scale of horizon crossing, in agreement 
with the WMAP5 central value. The normalization of the power spectrum at this point, 
^/V = 2 X 10~^, implies that the curvature scale is given by A; = 0.015 rrip. Eq. ( 4.15| ) shows 
that the radion mass is 



rrir 



Ve^k = 1.2 X 10"^ mp = 3 X 10^^ GeV. (4.35) 



From fig. H the scale of inflation mj = V^l'^ is given by = 3 x 10 ^ • dk'^nip in the flat 
region of the potential, hence 

mi = 2.7 X IQ-^mp. (4.36) 

The first slow roll parameter is e = 1.1 x 10~^, giving a small tensor-to-scalar ratio 
of r = 1.8 X 10^'^. This is related to the fact that the canonically normalized infiaton 
field does not change by an amount much greater than the Planck scale [^]; in our model 
A(j) = -^nipAip = Gnip. This is consistent with a refined version of the Lyth bound |37], 

A(j) = 6Mpri/^ (where Mp = VSTrnip). 

We investigated the effect of varying the parameters a and f3 away from the fiducial 
values a = 10~^ and /3 = —3. Holding /3 fixed and varying a, we find that for small a the 
spectral index asymptotes to its maximum value of Ug = 0.963 (figure |5|). Figure ^ shows 
that the curvature scale k, determined by the normalization of the power spectrum, follows 
a power-law relation k ~ a~^^'^ for very small values of a. Recall that the prefactor of the 
infiaton potential is proportional to k'^rrip] however the actual scale of infiation does not rise 
despite the increase in k; the e slow roll parameter asymptotes to 1.2 x 10~^ at the smallest 
values of a, which is nearly the same as its value at a = 10^'^. 

The corresponding results of varying /3 at fixed a = 10"'' are shown in figures ^ and ^. 
We see that it is possible to achieve a higher value of the spectral index than by varying only 
a; the maximum value is Ug = 0.99 at /3 = —3.1. 
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Figure 5: Spectral index versus logj^Q a for 
(3 ~ —3 model. 



1 




P 

Figure 7: Spectral index versus (3 for a — 
10~^ model. 




Figure 6: Log of curvature scale k versus 
logj^o a for /3 = —3 model. 




P 

Figure 8: Log of curvature scale k versus /3 
for a = 10~^ model. 



4.4 Model with larger tensors 

In the previous model, constructed from brane potentials consisting of just two exponentials, 
there was not enough freedom to find examples with a tensor-to-scalar ratio r exceeding 
10^^. However with the addition of just one more exponential term, we can overcome this 
limitation. Consider the Planck brane potential 

Vq = Aoe"«*/^5 _ Aoe'^o*/'^^ _ aQe^o^/^'^ (4.37) 
The back-reaction functions / and g are given by 

/ = -1 + e'^ (Ao e-"o ^ - ,5o e"* ^ + f^o e"^° + 5 (Aq oq 6""° ^ - <5o /3o e"^" + J^o <5o e"-^" ^) ) 
e"^ - Ao ao e""" ^ + Aq /3o e"'''' ^ - Aq e"^" ^ 



-Ao ao^e-"o Aq /Jo^e"''" - JIq Ao^e-^o ^ 



(4.38) 



which lead to the radion potential ( 4.11 ), proportional to e ^"^{f — g). 
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Solving the unperturbed junction conditions similarly to eq. ( p.18 - 2.19 ), we obtain 



An 



1 



Ao 



j3 — a 
1 



/3- 



a + 



{(3-6) 



{a -5) 



where we have reparametrized 5q = 5 + 1 similarly to ( 4.14 ). This model therefore has two 
parameters in addition to the previous one, 5 and Qq, the latter of which we find convenient 
to exchange for Qo/2/^|A; = ha. It will turn out that interesting values of Q,o/2fi'^k are of 
order a, so that h is of order unity. 

4.4.1 Analytic treatment 

We have seen that taking a ^ 1 is needed 
to get acceptable values of the spectral in- 
dex. The potential can be approximated by 
a simpler expression in this regime, that al- 
lows us to make analytic predictions for the 
tensor-to-scalar ratio r. By further tunings 
of parameters, we can achieve a potential 
that is nearly linear during inflation,^ lead- 
ing to a large e parameter and rj = 0. This 
is the best case for getting a large tensor 
signal from the model. Linear behavior in 
(/? occurs if some of the exponents in are 
quite small. Taking 5 is close to —3 realizes 

this possibility. 

First consider the small-a approximation to the potential, which can be written as 

l-6h- 3/ie-('^+3)'^ - e-3'^(l + - h{5 + 3)) 




Figure 9: Potential of extended model for h = 
0,1,2, and a = 1.5 X 10"^ ^ = -3, (5 = -2.9. 



Vr = 9k^mf,a 



3 + ae3^(4(l - 5h) - 3h{l + 6)^e-(^+^)v) 



(4.40) 



This is a good approximation in the inflationary region, between ip = and the maximum 
of the potential. To see the turnover of the potential at large ip, it is important to keep the 
term of 0{a) in the denominator; neglecting it leads to a potential that levels out at large (p 
instead of having a maximum. Next consider the case where 0<(5 + 3<Clso that the terms 
of order e~^'^ in the numerator are subdominant. This leads to 

1 - 5h - 3he-^^+^^^' 



Vr = 9k 



2 2 

nipa 



3 + aVsh{(p) 



(4.41) 



The term aVshif) is 0{1) near the maximum of the potential, but it becomes subdominant 
away from the maximum, and the 93-dependence is dominated by the numerator. In this 

^See ref. for another recent example leading to a linear potential, or ref. for other string-derived 
potentials supposed to be valid at Planckian field values. 
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region, since |5 + 3| <^ 1, the dependence is in fact linear. We can write ~ (1 — ^ + 8^99), 
where ^ = /i(3 + 6). This potential has slow-roll parameters t] = and 



I I m 



Inflation ends when e ^ 1, i.e., 



(1-^ + 3^(^)2 

2^-1 



3^ 



(4.42) 



(4.43) 



It is straightforward to solve the slow-roll equation of motion for the linear potential. 
Setting d-n/dN = in eq. ( |4.33| ), we have dip/dN = —{mp/fj,4)'^V^/Vr. This gives a simple 
quadratic function for N{(p). Using ( [4.43| ), one can find (f as a function the number of 
e-foldings before the end of inflation, 



(4.44) 



Notice that with A'^ = 60 and ^ = 1, 99(60) = 5, in agreement with the position of the 
maximum of the potential in fig. |^. Substituting this into (4.42) we get a simple result in 
which the dependence on the slope parameter ^ drops out completely: 



1 



4A^ + 1 



0.0044 



(4.45) 



where we have taken N = 56.5 to correspond to the COBE scale. This gives a tensor-to-scalar 
ratio of 



r = 16e = 0.07 

which is somewhat higher than the projected sensitivity of the Planck experiment, rn 
pH]. Moreover the spectral index is 



1 - 6e = 0.974 



(4.46) 



0.05 



(4.47) 



The prediction that r is independent of the parameters of the potential in this regime 
depends on the assumption that the linear behavior in the numerator of ( [4.41 ) dominates over 
the (/9-dependence from Vshif) in the denominator. In the next subsection we will quantify 
this by solving the equations of motion numerically using the full potential. 

Using eqs. ( [4.44 - 4.45 ), the power spectrum ( 4.34 ) in this model is P = 3{150tt'^)~^ a^k"^"^^ 
(AN + 1)^/2^ allowing us to determine ^/ak through the normalization: 



^/ak = 7.6 C"^/^ X 10"^ m. 



(4.48) 



The radion mass differs from that in ( [4.35 ) only by an extra factor of ^/T+^ ~ 1. 
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4.4.2 Numerical analysis 



Numerical analysis shows that the linear behavior of the potential is fully achieved for some- 
what smaller values of a ^ 10~^^ than we considered previously. For larger values of a, the 
linear shape is not fully realized and the tensor ratio does not reach its maximum value of 
r = 0.07. This is illustrated in figures ^ and 11, which show the dependence of r and 



n 



s on 

12 



h and 6 respectively, when a = 1.5 x 10 ^ . On the other hand, by decreasing a to 1.5 x 10 
and tuning 6 closer to —3, the predictions ( 4.46 ) and ( |4.47D based on the linear potential are 



borne out, both for r and for n^. In all cases, values of h{3 + 6) of order unity are needed to 



et the maximum tensor signal. 
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Figure 10: Tensor ratio and 1 — versus h 
for a = 1.5 X 10-^ /? = -3, (5 = -2.9. 




Figure 11: Tensor ratio and 1 — versus 8 
for a = 1.5 X 10"^ /3 = -3, /i = 3 
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Figure 12: Same as fig. 10, but for a = 1.5 x 
10-12 and 8 = -2.99. 



Figure 13: Same as fig. 11, but for a = 1.5 x 
10-12 and h = 100. 



It is encouraging that the model is able to produce values of r that exceed the minimum 
value of r = 0.05, which is estimated to be at the threshold for detection by the Planck 
satellite It is novel that we find such a regime in a hilltop potential [^], for which 

the tensor ratio is typically unobservably small, r < 0.002. We have checked that for the 
examples shown here, the maximum value of /, g in the inflationary region is / = 0.02, so the 
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condition f,g<^lis still satisfied and the 0{f'^,g'^) corrections to our approximation for the 
potential are under control. 

Consistent with the analytic prediction (|4.4^ ), the curvature scale k starts to exceed nip 
in the large-tensor models that have the smallest values of a. For example with 1.5 x 10^^^, 
6 = —2.99 and h = 200, we obtain r = 0.066 and k = lArrip. Examples with nearly as large 
tensors exist with k < nip, such as a = 5 x 10~^^, 6 = —2.985, h = 100, yielding r = 0.06 
and k = 0.8 m^. However it is the combination ^/ak that appears in the potential, and this 
remains subPlanckian, so the fact that k itself may exceed the Planck scale can be regarded 
as merely an artifact of our parametrization. 



5. Coupling of radion to standard model 

An interesting feature of this model that distinguishes it from most other models of inflation 
is that the couplings of the inflaton to the standard model are exactly specified, since the 
radion is a component of the higher-dimensional metric. This enables us to address the details 
of reheating in a way that is usually not possible without making additional assumptions. 
Moreover, we can compute the quantum corrections to the radion potential to check whether 
its exotic form is radiatively stable. 

The most natural situation is that the SM fields are localized on the TeV brane. As an 
example, consider a Higgs field H. Using the metric ( 3.10| ) and the conformal transformation 



(4.5), we find that in the Einstein frame, its coupling to is 

where we have renormalized the field and the mass to absorb warp factors so that mn is the 
already-warped mass, nin ^ TeV. For small ip, f{'p) ~ ay? and f2{'-p) ~ a/^v' which in the 
inflationary models we have considered is a negligible correction to the leading ip dependence 
in the exponents, giving the usual result that the radion couples to the trace of the stress- 
energy tensor. 



5.1 Reheating 

Using the conventional theory of perturbative reheating, one estimates the reheat tempera- 
ture as T.,. ~ (7* '^^^{TMpY/'^ where F is the decay rate of the inflaton and g^, the number of 



relativistic degrees of freedom. Using the interaction (5.1) and recalling the relation (4.8) be- 



tween and the canonically normalized radion (/> gives a rate of order F ~ g^^m^j / {IQirm^nir) , 
since there are ~ g.^ additional degrees of freedom besides the Higgs into which the radion 
can decay. For niH ~ 100 GeV, this leads to a reheat temperature that is too low even for 
nucleosynthesis, ~ 10~^ eV. 

However, there exists a much more efficient channel for reheating, by decay of the radion 
into gauge bosons. Although the radion does not couple to gauge bosons at tree level due to 
the tracelessness of their stress energy tensor, the conformal anomaly induces a coupling at 
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one loop, of the form 

^^^f^^^"" (5-2) 

where a is the fine structure constant. An analogous term is present for gluons and the 
electroweak vector bosons. The decay rate due to this operator is much faster than that 
due to ( ^.iD since the derivatives are of order mr/2 as opposed to mn-'^ The decay rate is 
approximately 

^ 32 2i3 7r3m2 ^ ^ 



Using the radion mass and a = 0.1 for QCD, this gives a reheat temperature of order 

Tr ~ 10^ GeV (5.4) 

which is high enough for electroweak baryogenesis, and low enough to avoid the gravitino 
problem. 

It is possible that reheating could be more efficient than indicated by this perturbative 
estimate, due to parametric resonance by the same coupling. Production of massless particles 
can be a particularly efficient form of preheating (see for example ref. [p]]). 

5.2 Radiative corrections to inflaton potential 

The radion potentials we have derived have exotic shapes compared to simple renormalizable 
potentials, and this has enabled us to produce distinctive signatures, like a larger tensor ratio 
than would be expected. One should always be concerned whether such peculiar potential 
shapes are radiatively stable. A second advantage of knowing the radion's couplings to matter 
is that we can address this question quantitatively. Consider the coupling { f>.l\ j to a Higgs 
field. To compute the contribution to the effective potential for ip, we can take <p to be 
constant and renormalize H to have a canonical kinetic term. Then the contribution to the 
Coleman- Weinberg potential is 

AV~.--^ (5.5) 

Since mjj- <C k'^nip (the prefactor of the radion potential), this is a negligible correction. Of 
course the more relevant contributions come from integrating out particles at some interme- 
diate scale mint, if there is new physics above the TeV scale. Even then however we are safe 
as long as mf^^ <^ k'^rrip. 

One correction that does not appear in the Coleman- Weinberg potential is the cosmolog- 
ical constant on the branes, a constant shift to the brane potentials ( 2.16| - ^T7 ). Such terms 



can be studied in the extended model with brane potential (4.37) in the limit where (^o = 



and by adding a similar term to the TeV brane. We have checked that our conclusions based 



of the potential with two exponentials in |4.3| are not very sensitive to the addition of the con- 



stant term. Fig. |lj shows the effect of the constant term on the potential: for h > 0, the flat 



*we thank Neil Barnaby for pointing out the possibihty of perturbative decay into gauge bosons 
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region of the hilltop is shortened, leading to a shorter period of inflation, while for h < —1 it is 
destabilized. The spectral index as a function of h is shown for the same potentials in fig. 15; 
values up to /i = 20 are compatible with the la WMAP5 allowed region, n<j = 0.963 it 0.015. 




Figure 14: Effect of constant addition 
to brane potential on inflaton potential, 
parametrized as eq. (4.17) with V,o/2fj,lk = 
ha, for a = 1.5 X 10~^, /3 = —3, and a range 
of values of /i = -10, 10, ... , 90. 




Figure 15: Spectral index as a function of h 
for the potentials of fig. 14. 



6. Discussion 

We have explored a new framework for inflation, where the radion in a warped 5D compact- 
iflcation is the inflaton. Our approach differs from the usual Randall-Sundrum setup in that 
we work in the linear dilaton background (a bulk scalar that varies linearly with the extra 
dimension) rather than AdSs. The bulk scalar stabilizes the radion by the Goldberger-Wise 
mechanism. Knowing its unperturbed solution exactly gives us greater computational control 
over the back-reaction when we perturb the radion away from the minimum of its potential. 
We derived an approximate analytic formula for the radion potential that depends mainly on 
the form of bulk scalar potential term on the Planck brane. 

Similar to the RS model, the hierarchy between the weak and Planck scales arises in 
our model without having to build it into the 5D Lagrangian. One can check that for the 
dimensionless parameter values we needed for inflation, the dimensionful brane potential 
parameters ( 2.18| , 2.19 ) are dominated by the terms Aq and Ai, both of order jj'^k for a warp 



factor consistent with kbo ~ 37 needed to solve the hierarchy problem. On the other hand, 
the mass of the radion does not get warped down to the weak scale; nir is flxed to be around 
10^^ GeV by the normalization of the CMB power spectrum. Moreover the couplings of the 
radion are not enhanced by the inverse warp factor; they remain Planck-suppressed. This 
tends to lead to a low reheat temperature, but also protects the inflaton potential from large 
radiative corrections. 

Choosing a sum of two exponentials as the brane potentials, we found that the radion 
potential suitable for inflation was generically of the hilltop form, but (with some tuning of 
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dimensionless couplings) having a shape that could be much flatter than the generic negative 
quadratic form. Expanding the ansatz on the Planck brane to a sum of three exponential 
enabled us to find potentials where the descent from the hilltop is linear and consistent with 
a large tensor contribution r = 0.07, detectable by Planck. 

An advantage of this model is that the couplings of the radion to matter on the branes can 
be computed explicitly. The radion couples to the trace of the stress energy tensor of a given 
particle, and to the one-loop conformal anomaly for states that are classically conformally 
invariant. We showed that this leads to reheating into gauge bosons by perturbative decay, 
with a reheat temperature of order ~ 10^ GeV. We did not explore the possibility of 
parametric resonance, which might significantly increase T^. 

A direction for further development is to see if our scenario can be realized in a string 
theory compactification. It is well-known that warped-throat-like linear dilaton solutions 
arise in the near horizon region of a stack of NS 5-branes [^] (see also section 14.1 of p3[). 
Moreover, in the same way as the Klebanov-Strassler throats are smoothly capped versions 
of anti-de Sitter throats near a stack of D3-branes, there exist known smoothly capped throat- 
like linear dilaton solutions |45], [^]. So far these throats and their compactifications have 
not been as widely studied as the analogous Klebanov-Strassler throat. 



In ref. [47| it was proposed that dynamically weakening gravity at early times could 
provide the low-entropy initial conditions that are needed for inflation to get started. In the 
Jordan frame our effective 4D action (4.3) has a time-dependent Newton's constant G{(j)) = 
(^g-</'//^4 when the radion is displaced from its equilibrium value, which could potentially 
realize this scenario. In this frame, 4D gravity is effectively weak when the radion is far from 
its stable minimum, which is also the condition needed for inflation. During inflation G{<j)) 
grows to its normal strength and then oscillates about this value as the radion decays. It 
would be interesting to further investigate the extent to which our model is compatible with 
this proposal. 

Another feature of our model is that the mass spectrum of general scalar perturbations, 
including the zero modes of the radion and bulk scalar field and the infinite tower of KK 
excitations, is completely known in terms of k and the parameters of the brane potentials. 
All modes can be normalized and their exact 4D effective actions determined. This presents an 
opportunity for constructing a model of assisted inflation [48| which we will explore elsewhere 
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A. Bulk field equations 

The ansatz of eq. (2.2) leads to the scalar field equation 



dt 



'N+3A+B 



d_ 

dy 



^N+3A~Bff^/ 



^N+3A+B 



d^ d<^> d^ ' 



and after substituting into the Einstein equations, Gmn = t^'^TpfN we have 



= 3e 



-2N 



+ Ab) - 3e~2^ {A" + 2A'^ - A B') 



G] = e"^^ { 2A + 3A^ + 2Ab - 2An - NB + B + B^ 



g-2B ^2A' + 3A'2 + N" + + 2AN' - 2AB' - N'B') 



Go5 



-2N 



A + 2A^- An) - 3e-2^ {A'^ + A'N') 



3e 



3 ( N'A + Ab - A A - A' 



(A.l) 

(A.2) 
(A.3) 

(A.4) 
(A.5) 



in the bulk where the indices are M,N = 0,1,2,3,5. The bulk contribution to the stress 
energy tensor comes from 



Tmn = 9mnV + dM^dN^ - ^ (d^^di^) qmn 
and we have the sources on the Planck and TeV branes: 



pM 



= e-^(°'*)5 (y) diag (^o, V^, Vb, 0) 
+ e-^(i'*)5 {y - 1) diag (Vi, Fi, Vi, Fi, 0) 



(A.6) 



(A.7) 



We do not consider the effects of adding matter to the branes, since we are primarily interested 
in the cosmology generated by the dynamical evolution of the radion coupled to the bulk. 

B. Kaluza-KIein excitations 



In this Appendix we provide details necessary for the computation of the KK mass spectrum 
summarized in Section [4.2| . We closely follow the methods of CGK [21| and KMP [^] and 
highlight differences in the zero mode solutions resulting from our different choice of stabilizing 
potentials on the branes. In particular, we have a second light mode due to not taking the 
stiff potential limit. 

Starting from the action of eq. ( |2.lD we expand to second order in the scalar fluctuations 



(^4.16|) to arrive at [|4| 



S = ^ I d^xv 



+ / d^xdy 



z"' 

-D+ — -- 



dy' 



z 



Jin 



2 3n 

n n * 



2nyl I 



(B.l) 
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where we define v and z as in eqs. ( 4.1^ . 4.19 ). The eigenvalue equations for separable solutions 
io V = Qj{x)vj{y) are in Schrodinger form 



-□ + 



z 
z 



(B.2) 



while the second term of (B.l) is not yet in a simplified form. To reach the compact form given 
by (|420| ) in |2| we use an equivalent set of eigenvalue equations involving F = Y2j Qji^)^jiy) 
which can be derived from the linearized Einstein equations. 

There are two nondynamical equations that we use to fix G and 6^ in terms of F. The 
first comes from recognizing that d^d,^ terms at linear order should vanish because there 
are no off-diagonal contributions in the stress energy tensor to source such terms. SG^i, = 
for these 4D components gives 

G = -2F (B.3) 



which has already been applied in The second constraint comes from the 5D linearized 

off-diagonal Einstein equation, (5G^5 = k^ST^^. Inserting ( p.3D it can be integrated to find 



(5$ 



l4 



%F' + 2n'F) 



(B.4) 



where a y-dependent "constant" of integration has been set to zero to ensure localization of 
the perturbations in x (the large dimensions). With these constraints, the system of linearized 
Einstein equations is then reduced to a single dynamical equation in the bulk found from the 



combination 



25GI 



F" + 3n'F' = DF 



(B.5) 



where one also eliminates V (eq. ( |2.7| )) using the background solutions and the property 
dV/d^ = 2fj,^^V of our bulk potential. Applying the definition ( 4.1S| ) and the constraint 
for eigenmodes OFj = —m'jFj, eq. (B.5) is seen to be equivalent to that derived from 



(I 

the action involving Vj (eq. ( p. 2D ). 

To determine the mass spectrum we must apply the linearized boundary conditions. At 
the Planck (-I-) and TeV (— ) branes. 



e ^ \y,-e 
"O \yi~e 



±2F^^^^V^{^i)\y^TS^^^fva<^^ 

T2Fv:i<^>,)y±6^v:'i^,)\y^ 



(B.6) 
(B.7) 



where primes on the potentials denote The ansatz for small fluctuations ( 4.16| ) makes 



( B.4 ) redundant with the first boundary condition. The second can be simplified using eq. 
(|3|) where we recall that the two would-be zero-modes get squared masses proportional to 



(Ti and o"2. If we also rewrite eq. (B^) in Schrodinger- like form 



2 T- 

-mjFj 



(B. 
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by defining 



so that 



F 



-3n/2jr 



(□ + m]) Qj = 
and if we moreover apply eqs. 



We have defined 



9 



3\y=yi 



dy^ 4' 

at the boundaries, then ( [B.?]) reduces to 

K{j k ho J~ j I y=y^ 



— — + - 



(B.9) 
(B.IO) 

(B.ll) 
(B.12) 



such that the index i identifies the brane position (y = yi) and j is an integer corresponding 
to the KK mode which we set to z = r, s when treating the light radion and bulk scalar 
modes. 

Solutions to ( [B.8D , which we provide below, and consequently also to ([B.2| ) when inserted 
back into the second order action (|B.1| ) , reduce to the effective actions of the KK modes given 
by eq. ( [4.2C ). We reproduce them here with the normalization coefficients given in terms of 



S = Y,Cj J d^x Qj [-□ - m| ] Qj 



Cj = /U5 / dy bo 



kbn 



2 

k 



yi=i 



yo=o 



(B.13) 
(B.14) 



We have also substituted for the bulk potential ( |2.7| ), eliminated the V^' terms in favour of 
perturbations using the boundary conditions ( [B.7| ), and applied the constraint equation ( |B.4D 
for (5$. 

The general solutions to the eigenvalue equations (B.IO) can be written similarly to eqs. 
for Fj 



Tj = Aj sin(\/Ajy) + Bj cos{^/Xjy) 



where as before 



2 

^3 = "^j - ^ 



4 



Az = — m: 

im?j > 9k'^ and <C k'^ respectively. 

Solving first for the KK modes, at y = yo = the constants are related by 

^ ^ _BjkKoj 



(B.15) 
(B.16) 

(B.17) 



(B.18) 
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At y = yi = 1 we get a transcendental equation which is solved to find the masses 



ky/Xj {kij — KQj) cos(y^ b) = [k'^Kij — Aj) sin( b) 



(B.19) 



The lightest KK excitation solution always gives Aq = with ttiq = The heavier modes 
must be found numerically except for the special cases of the stiff potential limit {\<Ji\ ^ 1), a 
massless radion (cjj = 0), fine-tuned brane potentials {ctq = cti) or that of our model where a 
large mass gap exists between the radion zero mode and the rest of the tower. In the former 
cases one solves sin(y^A^ 6) = to determine the masses. In our model when |cJo| <C |cri| ^ 1 
we find cos{^f\j 6) ~ and 



m 



4 



1\ TT 

^^2)kbQ 



n 2 



(B.20) 



We next turn our attention to the zero-mode solutions. Following the same reasoning 
usmg eq. ( pl^ ) and ( pill) at y = 2/0 = 0, we find the relation between the constants 



(B.21) 



At y = 2/1 = 1 we get the equation that determines m^, 



A^ + k^J \z{kqz - KxzYx - k KozKiz 



(B.22) 



where 



1 _)_ g-2v^fcbo 



1 



2v^fc6o 



(B.23) 



While eq. (B.22) does not yield a closed- form solution for the masses in general, one can find 



analytic results in some approximations. Expanding kyfXz — ^| for <C k"^ and 

ignoring the small exponential terms we can solve for the masses of the two light modes as 



ml ~ 



2(ToA;^ (71 fc^ 



1 



£0 

3 



1 + £1 
^ 3 



(B.24) 



in agreement with eq. ( 4.15| ). 



References 

[1] P. Horava and E. Witten, "Eleven-Dimensional Supergravity on a Manifold with Boundary," 
Nucl. Phys. B 475, 94 (1996) [arXiv:hep-th/9603142]. "Heterotic and type I string dynamics 
from eleven dimensions," Nucl. Phys. B 460, 506 (1996) [arXiv:hep-th/9510209]. 

[2] L. Randall and R. Sundrum, "A large mass hierarchy from a small extra dimension," Phys. Rev. 
Lett. 83, 3370 (1999) [arXiv:hep-ph/9905221]. 

[3] W. D. Goldberger and M. B. Wise, "Modulus stabilization with bulk fields," Phys. Rev. Lett. 
83, 4922 (1999) [arXiv:hep-ph/9907447]. 



- 32 - 



[4] N. Arkani-Hamed, S. Dimopoulos, N. Kaloper and J. March-Russell, "Rapid asymmetric 

inflation and early cosmology in theories with sub-millimeter dimensions," Nucl. Phys. B 567, 
189 (2000) [arXiv:hep-ph/9903224]. 

[5] N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, "The hierarchy problem and new dimensions 
at a miUimctcr," Phys. Lett. B 429. 263 (1998) [arXiv:hcp-ph/9803315]: I. Antoniadis, 
N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, "New dimensions at a millimeter to a Fermi 
and superstrings at a TeV," Phys. Lett. B 436, 257 (1998) [arXiv:hep-ph/9804398]. 

[6] J. M. Cline, "Inflation from extra dimensions," Phys. Rev. D 61, 023513 (2000) 
[arXiv:hep-ph/9904495] . 

[7] R. Sundrum and C. M. Wells, "Warped Hybrid Inflation," [arXiv:0909.3254 [hep-ph]]. 

[8] A. Vilenkin, "Topological inflation," Phys. Rev. Lett. 72, 3137 (1994) [arXiv:hep-th/9402085]. 

[9] A. D. Linde and D. A. Linde, "Topological defects as seeds for eternal inflation," Phys. Rev. D 
50, 2456 (1994) [arXiv:hep-th/9402115]. 

[10] R. C. Myers, "New Dimensions For Old Strings," Phys. Lett. B 199, 371 (1987). 

[11] E. Kiritsis, C. Kounnas and D. Lust, "A Large class of new gravitational and axionic 
backgrounds for four-dimensional superstrings," Int. J. Mod. Phys. A 9, 1361 (1994) 

[arXiv:hep-th/9308124]. 

[12] J. Polchinski, "String Theory, Vol. 1, An Introduction to the Bosonic String," Cambridge 
University Press (2000); sec cqs. (3.7.20) and (3.7.25). 

[13] O. DeWolfe, D. Z. Freedman, S. S. Gubser and A. Karch, "Modeling the fifth dimension with 
scalars and gravity," Phys. Rev. D 62, 046008 (2000) [arXiv:hep-th/9909134]. 

[14] P. Breitenlohner and D. Z. Freedman, "Stability In Gauged Extended Supergravity," Annals 
Phys. 144, 249 (1982); "Positive Energy In Anti-De Sitter Backgrounds And Gauged Extended 
Supergravity," Phys. Lett. B 115, 197 (1982). 

[15] P. Kanti, S. Lcc and K. Olive, "Stable, time dependent, exact solutions for brane models with 
a bulk scalar field". Phys. Rev. D 67 024037 (2003) [arXiv:hcp-ph/0209036]. 

[16] N. Kaloper, "Bent domain walls as brane worlds," Phys. Rev. D60, 123506 (1999). 

[arXiv:hep-th/9905210]. 

[17] J. M. Cline and H. Firouzjahi, "Brane world cosmology of modulus stabilization with a bulk 
scalar field," Phys.Rev.D 64 023505 (2001), [arXiv:hep-ph/0005235]. 

[18] D. Langlois and M Rodriguez-Martinez, "Brane cosmology with a bulk scalar field," Phys. Rev. 
D 64 123507 (2001), [arXiv:hep-th/0106245]. 

[19] K. Koyama and K. Takahashi, "Primordial fluctuations in bulk inflaton model," Phys. Rev. D 
67 103503 (2003). [arXiv:hep-th/0301165]. 

[20] K. Koyama and K. Takahashi, "Exactly solvable model for cosmological perturbations in 
dilatonic brane worlds," Phys. Rev. D 68 103512 (2003), [arXiv:hep-th/0307073]. 

[21] C. Csaki, M. L. Graesser and G. D. Kribs, "Radion dynamics and electroweak physics," Phys. 
Rev. D 63, 065002 (2001) [arXiv:hep-th/0008151]. 



- 33 - 



[22] J. M. Cline and H. Firouzjahi, "Five-dimensional warped cosmological solutions with radius 
stabilization by a bulk scalar," Phys. Lett. B 495 ,271 (2000) [arXiv:liep-th/0008185]. 

[23] C. Csaki, M. Graesser, L. Randall and J. Terning, "Cosmology of brane models with radion 
stabilization" Phys.Rev. D 62 045015 (2000). [arXiv:hep-ph/9911406]. 

[24] L. Kofman, J. Martin and M. Peloso, "Exact identification of the radion and its coupling to the 
observable sector," Phys. Rev. D 70, 085015 (2004) [arXiv:hep-ph/0401189]. 

[25] V. F. Mukhanov, "Gravitational instability of the universe filled with a scalar field", JETP 
Lett. 41, 493 (1985) [Pisma Zh. Eksp. Teor. Fiz. 41, 402 (1985)]. 

[26] M. Sasaki, "Large Scale Quantum Fluctuations In The Inflationary Universe," Prog. Theor. 
Phys. 76, 1036 (1986). 

[27] P. Brax, C. van de Bruck, A.C. Davis, and C.S. Rhodes, "Cosmological evolution of brane world 
moduU," Phys. Rev. D 67, 023512 (2003) [arXiv: hep-th/0209158]. 

[28] J. Lesgourgues and L. Sorbo, "Goldberger-Wise variations: Stabilizing brane models with a 
bulk scalar," Phys. Rev. D 69, 84010 (2004) [arXiv: hep-th/0310007]. 

[29] J. J. Blanco-Pillado et al, "Racetrack inflation," JHEP 0411, 063 (2004) 
[arXiv:hep-th/0406230]. 

[30] L. Boubekeur and D. H. Lyth, "Hilltop inflation," JCAP 0507, 010 (2005) 

[arXiv:hcp-ph/0502047] . 

[31] J. J. Blanco-Pillado et al, "Inflating in a better racetrack," JHEP 0609, 002 (2006) 
[arXiv:hep-th/0603129]. 

[32] Ph. Brax, S. C. Davis and M. Postma, "The Robustness of Ug < 0.95 in Racetrack Inflation," 
JCAP 0802, 020 (2008) [arXiv:0712.0535 [hep-th]]. 

[33] K. Kohri, C. M. Lin and D. H. Lyth, "More hilltop inflation models," JCAP 0712, 004 (2007) 
[arXiv:0707.3826 [hep-ph]]. 

[34] C. P. Burgess, J. M. Cline and M. Postma, "Axionic D3-D7 Inflation," JHEP 0903, 058 (2009) 
[arXiv:0811.1503 [hep-th]]. 

[35] D. H. Lyth and A. Riotto, "Particle physics models of inflation and the cosmological density 
perturbation," Phys. Rcpt. 314, 1 (1999) [arXiv:hep-ph/9807278]. 

[36] D. H. Lyth, "What would we learn by detecting a gravitational wave signal in the cosmic 

microwave background anisotropy?," Phys. Rev. Lett. 78, 1861 (1997) [arXiv:hcp-ph/9606387]. 

[37] G. Efstathiou and K. J. Mack, "The Lyth Bound Revisited," JCAP 0505, 008 (2005) 
[arXiv:astro-ph/0503360] . 

[38] L. McAllister, E. Silverstein and A. Westphal, "Gravity Waves and Linear Inflation from Axion 
Monodromy," [arXiv:0808.0706 [hep-th]]. 

[39] M. Cicoli, C. P. Burgess and F. Qucvcdo, "Fibre Inflation: Observable Gravity Waves from IIB 
String Compactiflcations," JCAP 0903, 013 (2009) [arXiv:0808.0691 [hep-th]]. 

[40] The Planck Blue Book, 

www.rssd.esa.int/SA/PLANCK/docs/Bluebook-ESA-SCI(2005)l.V2.pdf 



- 34 - 



[41] N. Barnaby, J. R. Bond, Z. Huang and L. Kofman, "Preheating After Modular Inflation," 
[arXiv:0909.0503 [hep-th]]. 

[42] C. G. . Callan, J. A. Harvey and A. Strominger, "Supersymmetric string solitons," 
[arXiv:hep-tli/9112030]. 

[43] J. Polchinski, "String Theory, Vol. 2, Superstring Theory and Beyond," Cambridge University 
Press (2000). 

[44] I. R. Klebanov and M. J. Strassler, "Supergravity and a confining gauge theory: Duality 
cascades and xSB-resolution of naked singularities," JHEP 0008, 052 (2000) 
[arXiv:hep-th/0007191]. 

[45] J. M. Maldacena and C. Nunez, "Towards the large N limit of pure N = 1 super Yang Mills," 
Phys. Rev. Lett. 86, 588 (2001) [arXiv:hep-th/0008001]. 

[46] A. H. Chamseddine and M. S. Volkov, "Non-Abelian solitons in N = 4 gauged supergravity and 
leading order string theory," Phys. Rev. D 57, 6242 (1998) [arXiv:hep-th/9711181]; 

"Non-Abelian BPS monopoles in N = 4 gauged supergravity," Phys. Rev. Lett. 79, 3343 (1997) 
[arXiv:hep-th/9707176]. 

[47] B. Greene, K. Hinterbichler, S. Judes and M. K. Parikh, "Smooth Initial Conditions from Weak 
Gravity," arXiv:0911.0693 [hep-th]. 

[48] A. R. Liddle, A. Mazumdar, and F. E. Schunck, "Assisted inflation," Phys. Rev. D58, 061301 
(1998), [arXiv:astro-ph/9804177]. 

[49] J. M. Cline and J. Trudeau, in preparation. 



- 35 - 



